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Abstract: The mass distribution of jets produced in hard processes at the LHC plays an 
important role in several jet substructure related studies involving both Standard Model 
and BSM physics, especially in the context of boosted heavy particle searches. We compute 
analytically the jet-mass distribution for both Z+jet and dijet processes, for QCD jets 
defined in the anti-A^ algorithm with an arbitrary radius R, to next-to-leading logarithmic 
accuracy and match our resummed calculation to full leading-order results. We note the 
important role played by initial state radiation (ISR) and non-global logarithms explicitly 
computed here for the first time for hadron collider observables, as well as the jet radius 
dependence of these effects. We also compare our results to standard Monte Carlo event 
generators and discuss directions for further studies and phenomenology. 
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1. Introduction 



Studies of jet substructure have become an area of great interest and much activity in 
the context of LHC phenomenology [1,2]. The primary reason for this relatively recent 
explosion of interest has been the observation that massive particles at the LHC can emerge 
with large boosts either due to the fact that one has an unprecedentedly large centre-of- 
mass energy y/s, which consequently enables access to a large range of particle transverse 
momenta, or due to the possibility that one may have very heavy new BSM particles that 
could decay to lighter Standard Model particles which thus emerge highly boosted. Such 
particles with energies or, equivalently, transverse momenta much greater than their mass 
Pt 3> m, would decay to products that are collimated and a significant fraction of the time 
would be recombined by jet algorithms into single jets. One is then faced with the problem 
of disentangling more interesting signal jets from those that are purely QCD background 
jets. In this context, detailed studies of jet substructure had been proposed some time 
ago [3] as having the potential to help separate signal jets (for instance, those arising from 
hadronic Higgs decays) from QCD background. More recent and detailed analyses [4-6] 
subsequently revealed the great potential offered by substructure studies in the specific 
case of boosted Higgs searches and paved the way for other substructure based studies in 
many different new physics contexts, which continue to emerge. 

However, one important issue that is a significant hinderance to the success of studies 
that try to exploit detailed knowledge of jet substructure is the complexity of the LHC 
as a physics environment, particularly in the context of strong interactions. For example, 
while looking for signal mass-peaks in the vicinity of a certain value of mass, m, one is 
inevitably led to consider the distribution of background QCD jets around that value. In 
fact a precise knowledge of jet-mass distributions and the related role of cuts on jet masses 
are extremely crucial for a wide variety of substructure studies. At the LHC one faces 
formidable obstacles as far as an accurate description of such spectra is concerned. Firstly, 
from the viewpoint of QCD perturbation theory, the most reliable tool at our disposal, one 
is confronted with large logarithms due to the multi-scale nature of the distribution at hand, 
which are as singular as — a™ln 2n_1 — , i.e double logarithmic, plus less singular terms. 
For boosted studies involving values pr 3> rnj such logarithms dominate the perturbative 
expansion and hence fixed-order tools like Nlo JET++ [7] , which are typically heavily relied 
upon for accurate QCD predictions, become at worst invalid or at best of severely limited 
use. To make matters worse, there are non-perturbative effects such as the underlying 
event and also the issue of pile-up, which leads to a contamination of jets resulting in an 
often very significant worsening of any perturbative description. 

In light of the above discussion, it is clear that the tools one can most readily use 
to estimate distributions of quantities such as jet masses are in fact Monte Carlo event 
generators. The parton showers encoded in these event generators derive from first princi- 
ples of QCD and offer a partial resummation of the large logarithms we mentioned before. 
They can be combined with fixed-order NLO results [8, 9] to yield descriptions that also 
describe accurately regions of phase space where the logarithms in question may not be 
entirely dominant. Moreover, they include hadronisation models as well as a model de- 
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pendent treatment of other effects such as the underlying event where the parameters of 
the model are extensively tuned to various data sets, to render them more accurate de- 
scriptions. While such tools are very general and hence of immense value in addressing the 
complexity of the LHC environment, the level of precision they offer may be considered a 
still open question. On the perturbative side the accuracy of the logarithmic resummation 
represented by parton showers is not clear. While leading logarithms (double logarithms in 
the present case) are understood to be guaranteed, at the level of next-to-leading or single 
logarithms (NLL level) the showers are not expected to provide a complete description. 
For example it is well known that parton showers work in the leading colour (large N c ) 
approximation as far as large-angle soft radiation is concerned, while the jet mass distri- 
butions we discuss here receive single logarithmic contributions from such effects which 
start already at leading order in a s and contain subleading colour terms. As we shall show 
single-logarithmic terms on the whole have a large impact on the final result and thus it is 
a little disconcerting to note that they will not be fully accommodated in current shower 
descriptions. 

Also the event generators we are considering have different models both for parton 
showers as well as for non-perturbative effects such as the underlying event and hadro- 
nisation and while it is always possible to tune the parameters in each event generator 
to experimental data, a comparison of the separate physics ingredients of each program 
often reveal large differences which, in certain cases, does not inspire much confidence in 
the accuracy of the final description, as far as QCD predictions are concerned. As one 
example of such differences we can refer the reader to the studies of non-perturbative ef- 
fects in jets at hadron colliders [10], where large differences were pointed out between the 
Herwig and Pythia underlying event estimates at Tevatron energies. A perhaps even 
more directly relevant example is the recent comparison by the ATLAS collaboration of 
their data on jet masses and shapes [11-13] to predictions from Herwig++ and Pythia 
which for the jet mass case, for example, do not agree very well with one another, with 
Pythia describing the data better. An understanding of the origin of such differences is 
certainly an important issue in order to gain confidence in the use of Monte Carlo tools for 
future LHC phenomenology. Hence the somewhat black-box nature of event generators is 
an issue which can makes sole reliance on these tools dangerous in the long run. 

Another avenue that can be explored in terms of theoretical predictions is the pos- 
sibility for analytical resummation of large logarithms. While the techniques needed to 
achieve such resummation currently apply to a more limited number of observables, while 
event generators are far more general purpose tools, where possible, resummation alleviates 
some of the present difficulties inherent in an event generator based approach. The typical 
accuracy of analytical resummation is usually at least NLL with some results having been 
obtained even up to the N 3 LL level [14]. It is also possible and straightforward in prin- 
ciple to match these resummed calculations to NLO estimates so as to have an accurate 
prediction over a wide range of observable values. Resummation has been a valuable tool 
in achieving high precision in several QCD studies ranging from pioneering studies of LEP 
event shapes [15] to current studies involving hadron collider event shapes [16]. In this 
paper we carry out such a resummation for the case of the jet mass in Z+jet and dijet 
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processes at the LHC. 

If one takes the study of jet-mass distributions as a particular case of studying shapes 
of jets produced in multi-jet events (see for instance Refs. [17, 18]) then it is clear that for 
substructure studies, resummation of the kind that we perform here for the jet mass, will be 
an important tool in achieving a precise description of the internal structure of jets. For the 
same reason, all the issues one encounters in the present study and the solutions we propose 
here to those problems, will also be of general relevance in the wider context of resummed 
calculations as substructure tools. There have in fact been several recent attempts to 
address the issue of resummation for jet masses and other substructure observables such as 
angularities [17-20]. In Refs. [17,18] it was proposed to study angularities of one or more 
jets produced in multijet events but calculations were only carried out for jets produced in 
e + e~ annihilation. The calculations carried out in these references also omitted important 
contributions at single-logarithmic level, the so called non-global logarithms [21,22], as 
was explained in some detail in Ref. [23]. Further calculations for jet-mass observable 
definitions in e + e" annihilation were carried out in Ref. [24,25]. References [19,20] on 
the other hand attempt to address the issue for the case of hadron collisions but employ 
calculations that are not complete to next-to-leading logarithmic accuracy, taking only 
into account the collinear branchings that generate the so called process independent jet 
function approximation to the resummed result. In this approach one does not treat the 
soft large-angle effects that arise from initial state radiation (ISR) or address the important 
issue of non-global logarithms or the dependence on the jet algorithm explained in Ref. [23], 
and hence cannot be considered sufficient for a reasonable phenomenological description of 
data. 

In our current paper we address both issues of ISR and non-global contributions and 
demonstrate their significance not just as formal NLL terms but also from the perspective 
of numerics and the accuracy of the final description for potential comparisons to data. 
We consider specifically the jet mass distribution of jets produced in two different hard 
processes: Z+jet production where it will be a background to the case of associated boosted 
Higgs production, with Higgs decay to bb and the case of jet production in dijet LHC events. 
We carry out a resummed calculation including the ISR contributions as a power series in jet 
radius R, while the non-global logarithms are calculated exactly at leading order (i.e. order 
a 2 s ) and then resummed in the leading N c approximation as was the case for DIS single 
hemisphere event shapes studied phenomenologically in Ref. [26]. We demonstrate that 
developing calculations for e + e~ variables and carrying them over to the LHC with neglect 
of process dependent ISR and non-global logarithms can yield large differences with the 
full resummation which correctly includes these effects. Moreover, our calculations have an 
advantage also over Monte Carlo event generators in that we retain the full colour structure 
of the ISR terms resorting to the leading N c approximation only for the non-global terms 
starting from order a^. The accuracy that we achieve in our resummed exponent should 
then be comparable to that which yielded a good description of DIS event shape data [26]. 
We also match our results to leading order QCD predictions so as to account for those 
terms which are not enhanced by large logarithms but may be important at larger values 
of jet mass, i.e. away from the peak of the distribution. The calculations of this paper are 
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valid for jets defined in the anti-A^ algorithm [27]. For jets defined in other algorithms such 
as the kt [28,29], and Cambridge- Aachen algorithms [30,31] the role of gluon self-clustering 
effects greatly complicates the single-logarithmic resummation (see e.g the discussions in 
Refs. [32-36]). For such algorithms analytical resummation to single logarithmic accuracy 
is still beyond the current state-of-the art and we shall hence not treat them. 

The calculations in the present paper also stop short of achieving the accuracy that 
was obtained for single-hemisphere DIS event shapes in one aspect. While we achieve 
the same accuracy as the DIS case for the resummed exponent, we do not yet obtain the 
NNLL accuracy in the expansion of the resummation as is achieved for most global event 
shapes in e + e~ and DIS in the leading N c limit [37] as well as in hadron collisions [16]. 
In other words our current resummation would not guarantee obtaining the a 2 s L 2 terms in 
the expansion, that arise from a cross-talk between a constant coefficient function a s C\, 
which corrects the resummation off just the Born configuration, and the a s L 2 term of the 
Sudakov form factors one obtains for jet masses. A proper treatment of such constant 
coefficients requires further work which we shall carry out in a subsequent article. At 
that stage we will also be in a position to carry out an NLO matching which, at least 
from the perturbative viewpoint, will give us an answer that will represent the state of the 
art for non-global observables. We do however estimate in this paper the possible effect 
of correcting for the coefficient function on our present predictions for the case of Z+jet 
production. In order to proceed to full NLL accuracy one would also need to understand 
non-global logarithms beyond the leading N c approximation but this is a much longer term 
goal. In the meantime we believe that the predictions we obtain here and certainly after 
forthcoming NLO matching will be sufficiently accurate so as to render them valuable for 
phenomenological studies. For the present moment we compare our resummed results to 
results from shower Monte Carlos and comment on the interesting features that emerge. 

We organise this article as follows: in the following section we outline the general 
framework for our resummed results indicating the separation between the global piece and 
the non-global terms. Following this, we derive in more detail the results for the global 
terms for both the case of Z+jet and dijet production. In section ||] we detail the results of 
our calculation for the non-global component at fixed-order and at all orders in the large 
N c limit. In section || we plot our final results for the case of Z+jet production, having 
carried out a leading-order matching and commented on the impact of various contributions 
to the resummed exponent such as ISR and non-global logarithms. We also discuss the 
expected effect on our results of a proper treatment of the coefficient C±, by treating its 
contribution in different approximations. Lastly, for the Z+jet case we compare our results 
to Monte Carlo estimates from a variety of event generators. In section |(| we discuss final 
results for the case of dijet production with matching to leading-order calculations. Finally 
in section |7| we arrive at our conclusions. Detailed calculations and explicit resummation 
formulae are collected in the Appendices. 

2. General framework 

The purpose of this section is to outline the overall structure of the resummed results that 
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we have computed for both jet production in association with a vector boson and dijet 
processes at hadron colliders. The notation we find most convenient to adopt is the one 
developed and used in Refs. [16,38] for the case of global event shape variables in hadron 
collisions. The extra ingredient specific to our calculation of jet masses is essentially that 
the observables we address are non-global, so that certain specifics shall of course differ 
from the case of global event shapes, which we shall highlight, where relevant. 

For the case of jet production in association with a Z boson we shall thus examine a 
distribution of the form 



where mj is the jet-mass squared of the highest px jet recoiling against the Z boson. 
For dijet production we shall instead adopt a different observable definition and study 
essentially the jet mass distribution averaged over the highest and second highest transverse 
momentum jets: 



where m 2 Jl and m 2 J2 are the squared masses of the highest and second highest pr jets, 
respectively. 

At Born level, for the Z+jet case, we have a single parton recoiling against the Z boson. 
If we restrict ourselves to a single-jet final state, then the jet-mass distribution for small jet 
masses is generated by soft and collinear emission off the hard Born configuration whose 
colour content is provided by the two incoming partons and the final state parton that 
initiates the jet. In contrast, for the case of dijet production soft and collinear emissions 
around the Born configuration, involve an ensemble of four coloured particles, with two 
incoming partons and two outgoing partons corresponding to the jets. 

While for global observables, such as event shapes, resummation off the hard Born 
configuration is all that matters, in the present case it is obvious that small jet masses can 
be produced in events with any jet multiplicity, that represent higher order corrections to 
the basic Born configurations we address. To restrict oneself to addressing just the Born 
configuration one can, for example, impose a veto scale pto, as suggested in Ref. [17]. How- 
ever, depending on the value of this scale one may then need to also resum the consequent 
logarithms involving the veto scale (see Refs. [17,18] for a discussion). Even if one chooses 
to adopt this procedure, the calculations we carry out and report in this paper shall still 
form the basis of the resummed answer, but will need to be modified to account addition- 
ally for the imposition of a veto. In the present article, we do not impose a veto but note 
that any additional production of non-soft, non-collinear particles (e.g the Z+2 jet correc- 
tion terms to the leading Z+jet process) will be associated with a suppression factor of 
o. s (pt) relative to the Born term, for each additional jet, where pt is the typical transverse 
momentum of the additional jet. This means that to the accuracy of our present paper 
(and indeed the accuracy of most current resummed calculations) we shall need to account 
for only the order a s correction to the Born term supported by a form factor involving 
only the double logarithmic (soft and collinear) component of the jet mass resummation. 
Thus, we never have to discuss, to our accuracy, the complex issue of soft wide-angle gluon 




(2.2) 
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resummation off an ensemble other than the Born configuration. The role of correction 
terms to the basic Born level resummation shall be discussed in more detail later in the 
article. 

Next, following the notation of Refs. [16,38] and denoting the Born kinematical con- 
figuration by B we write, for a fixed Born configuration, the cross-section for the squared 
jet mass to be below some value vp\, as 

d^ s \v) f 2 d 2 a^ 



/d a ( > 
d " ,2 '«^ e( ^- m 5>- <2 - 3) 



dB j „„„„ v 

The label 5 corresponds to the relevant production channel at Born level, i.e the flavour 
structure of the underlying 2—^2 Born process. We have also introduced the dimensionless 
variable v = rrtj/p^, with px the transverse momentum of the measured jet. One can 
integrate over the Born configuration with a set of kinematical cuts denoted by %{B) to 
obtain the integrated cross-section 

E«W = /*^«(5), (2.4) 

where, as should be clear from the notation, is the fully differential Born cross-section 
(i.e the leading order cross-section at fixed Born kinematics) for the subprocess labelled by 
5. We can then sum over Born channels 5 to obtain S(v) the integrated jet mass cross- 
section. However, the above result differs from the case of global observables considered in 
Refs. [16,38], in that it is correct only as far as the resummation of logarithms off the Born 
configuration is concerned and not at the level of constant terms which can arise from the 
higher jet topologies we mentioned before, that are not related to the Born configuration. 
When addressing the issue of constant corrections we shall thus need to account in addition 
to the above, for jet production beyond the Born level. For the present we focus on the 
basic resummation and hence work only with the Born level production channels as detailed 
above. 

Follwing Ref. [38] (for d <C 1) we then write 

The resummation is included in the function fjp and has the usual form [15]: 

ftp = ex P [Lgi(a 8 L) + g 2 {a s L) + a s g 3 (a s L) H ] (2.6) 

where gi, g 2 and g 3 are leading, next-to-leading and next-to-next to leading logarithmic 
functions with further subleading terms indicated by the ellipses and L = In ^ . 

For the observable we study here, namely non-global jet mass distributions, the func- 
tion gi is generated simply by the time-like soft and collinear branching of an outgoing 
parton and depends only on the colour Casimir operator of the parton initiating the jet, 
while being independent of the rest of the event. The function g 2 is much more complicated. 
It has a piece of pure hard-collinear origin, which, like the leading logarithmic function g±, 
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only depends on the colour charge of the parton initiating the jet and factorises from the 
rest of the event. In the collinear approximation, combining the soft-collinear terms of 
gi and the hard-collinear terms included in <?2 we recover essentially the jet functions first 
computed for quark jets in e + e~ annihilation in [15]. However, for complete single logarith- 
mic accuracy one has to consider also the role of soft wide-angle radiation. The function 
<72 receives a pure soft large-angle contribution also due to emissions from hard partons 
other than the one initiating the jet. For the Z+jet case this piece would be generated 
by coherent soft wide-angle emission from a three hard parton ensemble, consisting of the 
incoming partons and the outgoing hard parton (jet). For the case of dijet production, we 
have instead to consider an ensemble of four hard partons and the consequent soft wide- 
angle radiation has a non-trivial colour matrix structure [39], as for global hadronic dijet 
event shapes. 

Other than the above effects, which are all present for global event shapes and which 
are all generated by a single soft emission, the function 52 receives another kind of soft 
contribution, starting from the level of two soft gluons. Since we are looking in the interior 
of a jet rather than the whole of phase space, our observable is sensitive to soft gluons out- 
side the jet region emitting a much softer gluon into the jet. While for a global observable 
such a much softer emission would cancel against virtual corrections, in the present case 
it makes an essential contribution to the jet mass, triggering single logarithms in the jet 
mass distribution. These single logarithms (non-global logarithms) cannot be resummed 
by traditional methods which are based essentially on single gluon exponentiation. In fact 
a resummation of non-global terms, valid in the large N c limit, which corresponds to solv- 
ing a non-linear evolution equation [40], can be obtained by means of a dipole evolution 
code [21]. We carry out such a resummation in this article but do not attempt to address 
the issue of the subleading N c corrections which are as yet an unsolved problem. Since 
non-global logarithms are next-to-leading and we are in fact able to obtain the full colour 
structure for them up to order a 2 s , it is only single logarithmic terms starting at order 
where one needs to use the large N c approximation. One may thus expect that for 
phenomenological purposes an adequate description of non-global effects will be provided 
by our treatment here, as was the case for DIS event shape variables studied in Ref. [26]. 

In the next section we shall generate the entire result, except for the non-global terms, 
which we shall correct for in a subsequent section. The results of the next section correspond 
to the answer that would be obtained if the observable were a global observable and include 
process independent soft and hard-collinear terms alluded to above, as well as a process 
dependent soft wide-angle piece, which also depends on the jet radius R. This soft wide- 
angle piece, which starts at order a s , is calculated with full colour structure whereas one 
would expect that in Monte Carlo event generators only the leading N c terms are retained, 
thus implying higher accuracy of the results we obtain here. 

3. The eikonal approximation and resummation 

In the current section we shall consider the emission of a soft gluon by an ensemble of hard 
partons in the eikonal approximation. In this limit one can consider the radiation pattern 
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to be a sum over dipole emission terms [41]. Our strategy is to calculate the individual 
dipole contributions to the jet mass distribution and then to sum over dipoles to obtain 
results for both the Z+jet case as well as the dijet case. While the sum over dipoles shall 
generate both the soft-collinear and soft wide-angle terms we mentioned in the preceding 
section, we shall need to extend our answer to include also the relevant hard-collinear 
terms. Once this is done, the only remaining source of single logarithmic terms will be 
the non-global contribution to the single-logarithmic function g2, which we shall address 
in detail in the following section. 

We shall consider the most general situation that we need to address with all dipoles 
formed by two incoming partons and two outgoing hard partons. Clearly, for the Z+jet case 
one would have only a single hard parton in the final state, with the other parton replaced 
by a massive vector boson and hence for this case we will exclude the dipole contributions 
involving the recoiling jet. 

The squared matrix element for emission of a soft gluon k by a system of hard dipoles 
is described, in the eikonal approximation, as a sum over contributions from all possible 
colour dipoles: 

\M&\ 2 = \M B ,s\ 2 E C v W vW ' t 3 - 1 ) 

(ij)es 

where the sum runs over all distinct pairs (ij) of hard partons present in the flavour 
configuration 5, or equivalently, as stated before, over all dipoles. The quantity l-Mg^j 2 
is the squared matrix element for the Born level hard scattering, which in our case has to 
be computed for each separate partonic subprocess 5 contributing to the jet distribution 
and contains also the dependence on parton distribution functions. The contribution of 
each dipole Wij is weighted by the colour factor Cy, which we shall specify later, while the 
kinematic factor Wij(k) is explicitly given by the classical antenna function 

™ ~ (»")£•*) • (3 - 2) 

In the above equation a s is defined in the bremsstrahlung scheme [42], and its argument 
is the invariant quantity nf- = 2(pi ■ k)(pj ■ fe)/(pi • Pj), which is just the transverse 
momentum with respect to the dipole axis, in the dipole rest frame. We note that in 
the eikonal approximation, as is well known, the Born level production of hard partons 
in the relevant subprocess 5, factorises from the production of soft gluons described by 
the antenna functions W. The squared matrix element |-Mg l( s| 2 essentially produces the 

quantity —A- while the W functions start to build up the exponential resummation factor 
fjp referred to in Eq. (j2.5|). In what follows below we shall focus on the various components 
of the resummation in more detail and in particular carry out the calculations from the 
individual dipole terms. 

3.1 Exponentiation: the Z+jet case 

We have mentioned above the antenna structure of soft gluon emission from a system of 
hard emitting dipoles. It is well understood by now that if one ignores configurations cor- 
responding to non-global logarithms (in other words those that stem from emission regions 
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where gluons are ordered in energy but not in angle), then to single-logarithmic accu- 
racy there is an exponentiation of the one-gluon emission terms described in the preceding 
section, as well as the corresponding virtual corrections that have the same colour and 
dynamical structure but contribute with an opposite sign so as to cancel the divergences 
of real emission. We note that for the case of Z+jet we are dealing with a hard parton 
ensemble with three partons, two incoming and one corresponding to the triggered jet. In 
this case the colour factors Cij = — 2(Tj.Tj) (with the Tj being SU(3) generators), that 
accompany the dipole contributions, can be straightforwardly expressed in terms of quark 
and gluon colour charges. Taking account of virtual corrections and the role of multiple 
emissions generating the jet mass one can write the result in a form that is familiar from 
the earliest studies of jet masses in e + e~ annihilation [15] 

AS) _ exp[-ft g - Tg^j] ,„ 
/s.giobai- r(l + K' s ) ' { ' 

where the subscript above denotes that we are considering the global term only i.e ignoring 
all non-global corrections to fg\ 

The function {—TV) represents the exponentiation of the single-gluon contribution after 
cancellation of real-virtual divergences, while 71' is the logarithmic derivative of 7Z, OlTZ 
to be evaluated to our accuracy simply by accounting for the leading logarithmic terms in 
71. The terms involving TZ' arise due to the fact that direct exponentiation only occurs 
for the Mellin conjugate of the variable v. To single-logarithmic accuracy one can invert 
the Mellin transform analytically by performing a Taylor expansion of the Mellin space 
result to first order and integrating over the Mellin variable, resulting in the form written 
above [15]. 

One has for 7Z$ the result: 

7l 5 = Y, [ °ij dk th dr, ^ Wij (k) e (v(k) - v) , (3.4) 
(ij)es 

where we have introduced the integral over the momentum of the emitted gluon k and the 
step function accounts for the fact that real-virtual cancellations occur below a value v 
of the normalised squared jet mass, while uncancelled virtual corrections remain above v. 
The function v(k) is just the dependence of jet-mass on the emission k. Letting the hard 
initiating parton have rapidity y and transverse momentum pt and denoting by kt, T] and 
<j) the transverse momentum rapidity and azimuth of the soft gluon k we have (when the 
hard parton and gluon are recombined to form a massive jet) 



v(k) = = ^ [cosh {r, - y) - cos 4,] + O (| ) (3.5) 

where we neglect terms quadratic in the small quantity k t . 

It now remains to carry out the dipole calculations for the Z+jet case. We have a 
hard process with two coloured fermions and a gluon or a three-hard-particle antenna, 
irrespective of the Born channel 5. Let us call Si the Born subprocess with an incoming 
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quark (or anti-quark) and an incoming gluon, which results in a final state coloured quark 
or antiquark recoiling against the Z boson. Labelling the incoming partons as 1 (fermion) 
and 2 (gluon) and the measured jet as 3 we have the following colour factors: 

C12 = N c , C 2 s = N c , C 13 = --L (3.6) 

For the remaining Born channel with an incoming qq pair we obtain the same set of colour 
factors as above but with an interchange of 2 and 3, to correspond to the fact that is always 
the quark-(anti)quark dipole which is colour suppressed. 

The calculation of individual dipole terms is carried out in Appendix [A]. We use the 
results obtained there to construct the final answer. Let us focus on the Born channel 5± 
corresponding to an incoming gluon and quark with a measured quark jet. In order to 
combine the various dipoles that contribute to the resummed exponent, we combine the 
pieces TZ{j computed in the appendix weighting them appropriately by colour factors: 

n Sl = C12 n l2 + Ci 3 (n\f + raff) + c 23 (n s 2 f + n c 2 f) . (3.7) 

Using the appropriate dipole results generated by using the eikonal approximation, 
from Appendix [A] and using the colour factors mentioned in Eq. (fOj) we obtain 

AT 2 — 1 /V 2 — 1 

n Sl = ^^n con - + N c n\f + -^-—n^ , (3.8) 

where we used the fact that = K c 2 f- = K co11 and K\f = K s 2 f = K soh . Writing 

the result in terms of the colour factors Cf and Ca and the explicit results for the various 
dipoles results in the following simple form: 



n Sl (v) = 2C F J 



a s (h,j) dk jj 
2vr kl 



a s {h,j) dkf j 





+ 2C F / a ; ^ In 

J 2vr kfj 

where kt j is the transverse momentum of the emitted gluon with respect to the jet. 

The above result represents the decomposition of the resummed exponent into a 
collinear piece contained in the first two lines of the above equation and a soft wide- 
angle piece. We have included in the collinear piece a term e -3 / 2 in the argument of the 
logarithm, which corrects the eikonal approximation for hard collinear splittings of the 
final state quark jet. This correction term emerges from replacing the IR singular (pole 
part) of the q — >■ qg splitting function, treated by the eikonal approximation, by the full 
splitting function. As one would expect this collinear piece, which also contains the lead- 
ing double logarithms, involves only the colour charge Cf of the parton that initiates the 
measured massive jet, in this case a quark. The remaining part of the result above is a 



- 11 - 



process dependent soft large angle piece that has a power series expansion in jet radius R, 
which we truncated at the R 4 term. We note that the R term emerges with a numerically 
small coefficient and shall make a negligible impact on our final results which can thus be 
essentially obtained by considering the 0(R 2 ) corrections alone. We also note that the 
calculation of the soft large-angle component of the result should mean that our results 
are more accurate than those obtained from MC event generators, which would only treat 
such pieces in a leading N c approximation. 

The case of the other subprocess, which generates a gluon jet in the final state, is 
totally analogous. The result for the resummed exponent is 

,» («*-*™*) e ( k \ - v ) e (# - Bf) 
27r Kj \ kt > J \Pt \ Pt J 



+ R 2 {^ F -^j' d ^ ( ^ + 0{R% (3.10) 

In order to achieve NLL accuracy, the remaining integrals in Eq. fl3.9|) and Eq. ( [3.10 ) must 
be performed with the two-loop expression for the running coupling. We obtain: 

K Sl = -C F (Lfi + f 2 + / coll , g ) - i? 2 /i.a. (c A + ^) + 0(R A ), 



■Rs,- -(' i ( LU + h + ./coil,,) - R 2 fi. a . [2C F -^)+ 0{R^ 



(3.11) 



with L = In ^- . Explicit results for the functions fi are collected in Appendix [C]. 

3.2 Exponentiation: the dijet case 

We now turn our attention to the process 

p(P 1 )+p(P 2 ) -> J( P3 ) + J(p 4 ) + X, (3.12) 

where we want to measure the mass of the two leading jets. For convenience, we fix the 
kinematics of the two (back-to-back, in the eikonal limit) leading jets, i.e. their transverse 
momentum p F and their rapidity separation |Ay|. The calculation of the dipoles in the 
eikonal limit proceeds in the same way as in the Z+jet case that we have previously 
analysed. The main difference is the more complicated colour algebra, that leads to a 
matrix-structure of the resummed result. The formalism to perform the resummation in 
the presence of more than three hard partons was developed in [39]. For each partonic 
sub-process we need to fix a colour basis and find the corresponding representations of the 
colour factors Tj.Tj. For each partonic subprocess, the resummed exponent takes the form 



AS) _ 1 \- 
^global" tlH 2^ tr 
J=3,4 



Hse -(Ss,j+y E g' s>j y s StJ e~ Ss ' J ^ B0 '^ + (Ay <-> -Ay) 

r (i + 205, j 



(3.13) 
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da {&) 

The matrices H$ correspond to the different Born subprocesses and txH = Jg . We note 
that the resummed expression in Eq. ( |3.13| ) is written in terms of exponentials that describe 
the colour evolution of the amplitude 1 , rather than of the cross-section as in the Z+jet 
case. We obtain 

Gs,s = (Lfx + f 2 + /coii, 3 ) + T 1 .T 2 / 1 . a .(27ri + R 2 ) 

+R 2 fi. a . QT 3 .T 4 tanh 2 ^ + i(Ti.T 3 + T 2 .T 3 ) 

1 e Ay 1 e~ Ay \ 4 

+ 2 Tl - T4 l + coshAy + 2 T2 - T4 l + coshAy ) + ° {R } ' 

Q &A = (Lfi + f 2 + /coii,4) + Ti.T 2 / La .(27ri + R 2 ) 



(3.14) 



+^/i.a. ^-T 3 .T 4 tanh^ + -(T1.T4 + T 2 .T 4 ) 

1 e Ay 1 e~ Ay \ 
+-T 2 .T 3 — + -T1.T3 + 0(R 4 

2 d l + coshAy 2 1 + coshAyy V 



where the functions fa are reported in Appendix |C] and, as before, L = hi-^-, Q' = diG- 
The collinear part of the result is diagonal in colour space, with a coefficient which is the 
Casimir of the jet. Large-angle radiation is instead characterised by a more complicated 
colour structure. We also note the presence of the imaginary phase due to Coulomb gluon 
exchange. We choose to work in the set of orthonormal bases specified in [43], to which we 
remind for the explicit expressions. As a result, all the colour matrices are symmetric and 
the soft matrix appearing in Eq. ( 3.13| ) is the identity Ss,j = 1. 



As an example, we report explicit results for the scattering of quarks with different 
flavours q(i)q'(j) — >• q(k)q'(l). We work in a normalised singlet-octet basis: 

N c 

In the t channel (Ay > 0), the hard scattering matrix is given by 



= Ml 4* ' (3 ' 16) 



In the literature this resummed expression is written in terms of an anomalous dimension F, where 
Q = T£ and £ is the appropriate evolution variable. 
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We have that T§ 



rp2 

A 4 



Ci? and the other colour matrices are 



Ti.T 2 = T 3 .T 4 



T1.T3 = T 2 .T 4 



T 2 .T 3 = Ti.T 




(3.17) 



3.3 The constant term C\ 

In order to achieve the NNLL accuracy in the perturbative expansion, that is common 
for most resummation studies of event shape variables, one must consider also the 0(a s ) 
corrections, which are not logarithmically enhanced in the small jet-mass limit and their 
cross-talk with double logarithmic terms arising from the Sudakov form factors. The con- 
stant terms can be expressed as : 

' "da® 



a s C[ 



(S) 



lim 

v-*o L 



J NLO 



(V) 



,(-5) 



lim 



(<5) 
NLO 



' ^NLL,a s ( V ) 

Umax ^ a (S) 



lim 

u->0 



dv 



-dv 



,(•5) 

J NLL,a 



(V) 



a NLO + lim 



dv 

da^ 
dv 



dv — S 



(S) 



NLL 



dv 



J NLL 



(3.18) 



If 5 is a partonic channel that is also present at Born level, then we can recover the usual 
definition of the constant term: 



(8) _ Ci 



(3.19) 



The general kinematic, colour and flavour structure of can be rather complicated. 
However, as discussed for global event-shapes in Ref. [16,38], one can simply multiply 
together this constant and the appropriate resummed exponent previously discussed, 
essentially because the only relevant terms at NLL originate by the product of cf^ times 
the double logarithms (soft and collinear) in the exponent, which do not depend on the 
colour flow in the hard scattering nor on the parton distribution functions. This is also true 
in the case of the jet-mass we are considering in this paper, with the further complication 
that we need to specify the flavour of the jet we are measuring, because quark and gluon jets 
will receive different suppression factors. In particular, new channels open up at relative 
order a s which are not related to the Born channels and hence obtaining the contribution 
of the constant terms separately for these channels is an involved exercise. While we leave 
the complete determination of C\ in different experimental set-ups, as well as an analysis 
of the dijet case, for future work, for the case of the jet mass in Z+jet, where we measure 
the mass of the hardest jet, we shall provide later an estimate of the contribution of the 
constant terms C\ to the resummed distribution. 
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4. Non-global logarithms 
4.1 Fixed order 

As we emphasised before, the jet mass that we study here, is a non-global observable which 
means that the results presented in the previous section are not sufficient to obtain the 
next-to-leading logarithmic accuracy that is common for event shape variables in e + e~ 
annihilation. One is also required to correct the results for the effect of soft wide-angle 
emissions arising from gluon branching; in other words correlated gluon emission as opposed 
to independent emission of soft gluons by the hard particle ensemble. Calculations involving 
correlated gluon emission have been carried out at fixed-order and all-orders (in the leading 
iV c limit), in the simpler cases of non-global event shapes in e + e~ annihilation (such as 
the light jet-mass) and DIS. Till date a full calculation even at fixed-order has not been 
carried out for hadron collisions, which we shall do below, in the context of the jet-mass 
distribution. We also note that in our previous work we carried out a calculation of non- 
global logarithms for jet masses of jets produced in e + e~ annihilation, in the limit of small 
jet radius R <C 1, which we argued could serve as a model for the calculation in the hadron 
collision case. However, in the current work we shall lift the approximation of small R 
meaning that our calculations should be useful for jets of any radius. 

To address the correlated emission term at leading-order we need to consider the case 
where one has two energy-ordered soft gluons k\ and k 2 such that, for instance, u\ S> u>2, 
where ui\ and U2 are the respective energies. In previous work involving e + e~ annihilation 
we have addressed this situation by using the fact that the emission of such a two-gluon 
system off a hard qq dipole is described by an independent emission term with colour factor 
Cp and a correlated emission term with colour factor CfCa- However, in the present case 
of multiple hard partons, the emission of a two-gluon system is in principle more involved, 
since there are several emitting dipoles. In practice, the structure of two-parton emission 
for the cases we study in this paper (i.e up to n = 4 hard legs) is known to be remarkably 
simple [44,45]. 

As an illustrative example we consider again the Z+jet case, with three hard legs 
and the Born channel S±, with an outgoing quark jet. The squared matrix element once 
again contains an independent emission piece, which contributes to the exponentiation of 
the single gluon result, described by our function 1Z. This leaves the correlated parton 
emission term which has the structure 

(2) 

where is the correlated two-gluon emission by a dipole ij, and which is the same as 
for the qq case studied in e + e~ annihilation. Hence the dipole emission and associated 
colour structure for a correlated two-parton system is precisely the same as for a single 
gluon emission [45]. 

Next we note that, as described for example in Ref. [45], a piece of the correlated 
two-parton emission contribution actually goes to build up the running coupling we have 
already considered in the exponentiated single gluon contribution 1Z. This piece comprises 
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of gluon splittings into equally hard gluons or into a qq pair, which together produce the 
leading term of the QCD /3 function 2 . This leaves us to consider only the soft part of the 
correlated emission S [46], which describes the production of an energy-ordered two-gluon 
system. For a global observable as is well known, this term produces no effect at single- 
logarithmic accuracy whereas in the present case it provides us with the first term of the 
non-global contribution. For a general dipole ij we can explicitly write 

Sij = 2C A T i .T j (A ab + A ba - A a A b ) , (4.2) 

where 

Aab = M Aa = fa") (4.3) 

(ik a )(k a kb)(hj) ' a (ik a )(jk a ) 

We note that this term is free of collinear singularities along the hard legs a and b due to 
cancellation between the various terms of Eq. ( f4.2|) . The remaining collinear singularity 
from the 1/ (fcj.^) , involving the soft gluons, shall turn out to be integrable for the non- 
global configurations considered below. Thus to obtain the leading non-global term for the 
jet mass under consideration, it suffices to carry out dipole calculations using the Eqs. ( fO^ ) 
and (|4.3| ), for each of the hard emitting dipoles, just as we did for the single gluon emission 
case. The results we obtain are mentioned and commented on below with the details on 
the calculation in Appendix ^. We report below the coefficients hj of the non-global single 
logarithms: 

J [dk^d^SijQk^jQk^j = C A Ti.T 3 I i3 In 2 ~. (4.4) 

• dipoles involving the measured jet 

As was demonstrated in our previous work [23], for such dipoles one obtains a re- 
markably simple result independent of R and dipole geometry: 

r 2 
3 



7T 

13 = hz = Ha = ^r- (4.5) 



This result is precisely the same as was obtained for hemisphere jet masses in e + e~ 
annihilation and implies that the result obtained does not depend on the position of 
the jet boundary, a feature that has also been observed in [47]. 

in-in dipole 

Here we are considering a situation where the emitting dipole legs are away from the 
interior of the jet region into which the softest gluon &2 is emitted. This situation is 
reminiscent of the much studied case of non-global logarithms in energy flow into gaps 
between jets [22,40]. In Ref. [22] for instance, energy flow into an inter-jet region (i.e 
a region between but away from the hard legs of an emitting dipole) was considered 
for specific choices of the geometry of this region, such as a rapidity slice and a patch 
in rapidity and azimuth. In the present case of course we are studying the same 



2 For a non-global observable such hard splittings will also give rise to non-global logarithms below the 
single logarithmic level, which are subleading from our point of view [25,34]. 
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problem but the region being considered here is a circle (77 — y) 2 + (f> 2 , centred on the 
measured jet. The result we obtain for the non-global contribution due to the in-in 
dipole was computed numerically but has the small R behaviour 



J12 «4[l.l7i? 2 - R 2 \n2R + 0{R i )] . (4.6) 

We note that the R 2 behaviour arises simply from integrating the emitted softest 
gluon &2 over the interior of the jet region while the In R behaviour is a reflection of 
the collinear singularity we mentioned above, between the two soft gluons k\ and &2- 
In practice the above small R approximation is not a good approximation to the full 
result obtained numerically, for larger values of R ~ 1, hence eventually we use the 
numerical answer rather than the form above. 

• in-recoil dipole 

Here again the situation is similar to the case of the interjet energy flow with the only 
difference from the previous dipole being the geometry of the hard dipole legs, which 
are now formed by an incoming parton and an outgoing parton at a finite rapidity 
with respect tothe beam. The result reads 

j 12 „ l i + 6 ' „ (1.17R 2 - R 2 In 2R) + 1 + - { n(Ay)R 2 . (4.7) 
(l + coshAy) 2V ; 1 + coshAy v yj y ' 

where k depends on rapidity difference Ay = y — y r and is evaluated numerically. In 
the limit Ay — > 00, n(Ay) — > and one recovers the previous case. The dipole I24 is 
easily obtained as 

hi = Iu(-Ay). (4.8) 

4.2 All-order treatment 

In order to perform a complete (NLL) resummation of non-global logarithms, one would 
need to consider the emission of a soft gluon from an ensemble of any number of gluons. 
This problem can be treated only in the large- N c limit [21, 40] 3 . For this study, we have 
adapted the dipole-evolution code used in [21] to perform the resummation in the case of 
jet-masses. 

The code developed in Ref. [21] handled the case of evolution off a hard primary dipole 
in the leading N c limit. The result for the non-global contribution S(t) was obtained by 
dividing the resummed result by the contribution of primary emissions alone, off the hard 
emitting dipole. For our work here we have a situation with several possible emitting 
dipoles. In this situation one has to resort to the large iV c limit in which one can treat the 
problem as independent evolution of only the leading colour- connected dipoles in the hard 
process. Detailed formulae can be found in Ref. [40]. 



3 An alternative approach that can be found in the literature consists in an expansion in the number of 
out-of-jet (in this case), gluons, keeping the full colour structure, see for instance [43,48-50]. 
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In the Z+jet case, which is simpler, we have noted that there is no visible difference 
arising from considering just the leading colour-connected dipoles in the hard process rel- 
ative to the case where one evolves all hard dipoles using the evolution code. We choose 
the latter option here hence write the full resummed expression, with the exception of the 
contributions coming from the constant terms C\, as 



= E / 



da (8) s 
rg /s,global/s,uon-global^H^)- 



(4.9) 



The resummation of non-global logarithms, including contributions from the colour sup- 
pressed hard dipoles, is encoded in the two terms 



Ah) 

J non- global 



Ah) 

J non— global 



exp 



exp 



7T" 



-CaCf— /l3(*) 



7T 



C A 



ci^-Mt) -C A [C F -^\ h 2 (R)fi2(t) 



(4.10) 
(4.11) 



where we have used that J13 = I23 = I34 = the other contribution I12 does depend 
on R and vanishes in the R — > limit, where one recovers the picture of jets evolving 
independently. As we stated before, taking the large N c limit of the non-global contribu- 
tions would amount to switching off the contribution from the colour suppressed dipoles 
or, equivalently, choosing Cf = Ca/2 in the above. This produces no significant difference 
to our results but the result written above has the advantage that it includes the the full 
contribution to 0(a%) non-global coefficient. We have defined 



fij(t) 



1 + {flijtf 



■if 



t 



1 



In l-2a s (pr)A)ln 



R- 



(4.12) 



The coefficients a^- , h l3 , 



C{j 3X6 obtained by fitting the functional form above to the nu- 
merical results from the large- iV c dipole-evolution code. Numerical results are reported in 
Appendix |^. 

Following a similar method, we obtain the corresponding result for the dijet case: 



E 

5,J=3,4 



dBtr 



-(Ss,j+7bS' Si j+S S) j) -Qs : j-1eQ' Si j-Ss : j 



+ (Ay <-> -Ay) 



r i + 25' 



5,J 



(4.13) 



with 



T§/ 13 (i) - Ti.T 2 Ji 2 CR)/i2(t) - T^TtluiR, Ay)f u (t) 



Ca/tS 
2 V 3 

-T 2 .T A I 2i (R,Ay)f 2i (t) 

^(yTl/ 13 (t) - Ti.T 2 ii 2 (12)/i2(t) - T 2 .T 3 I U (R, Ay)f 14 (t) 
-Ti.T 3 /24(12 ) A I /)/ 2 4(*) 



(4.14) 

As before the above expressions capture the full colour structure of the non-global contri- 
bution at 0(0?), but beyond that are valid only in the large- N c limit. 
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5. Z+jet at the LHC 



In this section we investigate the numerical impact of the different contributions which are 
relevant in order to achieve NLL accuracy. We decide to study the differential distribution 
where £ = y/v = , so as to study the jet mass distribution directly rather than 
the squared jet mass. We also find it useful to work with a dimensionless ratio to better 
separate soft physics contributions. In fact, a fairly large value of the jet mass can be 
generated by the emission of a very soft gluon, if the transverse momentum of the hard jet 
is large, while small values of £ always correspond to the emission of soft and/or collinear 
gluons. If not stated otherwise, we normalise our curves to the Born cross-section. We use 
the matrix element generator Comix [51], included in Sherpa [52], to produce all the tree- 
level cross-sections and distributions. We consider proton-proton collisions at y/s = 7 TeV 
and we select events requiring px > 200 GeV; the Z boson is produced on-shell and does 
not decay. Jets are defined using the anti-fet algorithm [27]. In our calculation we use the 
set of parton distribution function CTEQ6m [53], with renormalisation and factorisation 
scales fixed at /ur = fj,p = 200 GeV, to ease the comparison with different Monte Carlo 



parton showers, which we are going to perform in Section 5.4 



5.1 Different approximations to the resummed exponent 

We start by considering different approximations to the resummed exponent /» ■ We 
present our results for two different jet-radii: R = 0.6, in Fig. [l], and R = 1.0, in Fig. |2| The 
blue curve corresponds to the most simple approximation to the NLL resummed exponent: 
the jet-function. This approximation correctly resums soft and collinear radiation as well 
as hard collinear, but does not capture all soft radiation at large angle. In particular this 
corresponds to neglecting terms that are suppressed by powers of the jet radius in the 
resummed exponent Eq. ( |3.11 ). These terms are included in the resummation of all global 



contributions (green line). We have checked that inclusion of 0(R 2 ) is enough because 
the 0(R 4 ) corrections are below the percent level even for R = 1.0. We stress that up 
to this point no approximation on the colour structure has been made, although we have 
checked that sub-leading colour corrections are small, once the collinear part has been 
properly treated. Finally, in the red curve we also take into account the resummation of 
non-global logarithms as described by Eq. ( |4.9| ). The first, 0(a 2 ), coefficient on the non- 
global contribution is computed exactly, while the subsequent resummation is performed 
using a numerical dipole-evolution code in the large- iV c limit. We note that the inclusion 
of 0(R 2 ) terms in the resummed exponent as well as non-global logarithms, noticeably 
corrects the simple jet-function picture, based on collinear evolution. The peak height is 
reduced by more than 30% for R = 0.6 and it is nearly halved for R = 1.0. The effect 
of non-global logarithms is reduced in the latter case, but the 0(R 2 ) corrections to the 
jet-function approximation become bigger. 
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Z+jet, R=0.6, p TJ > 200 GeV 
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Figure 1: Comparison between different approximations to the resummed exponent: jet functions 
(blue), with full resummation of the global contribution (green) and with non-global logarithms as 
well (red). The jet radius is R = 0.6. 




Figure 2: Comparison between different approximations to the resummed exponent: jet functions 
(blue), with full resummation of the global contribution (green) and with non-global logarithms as 
well (red). The jet radius is R = 1.0. 



5.2 Matching to fixed order 

We now turn our attention to obtaining a jet-mass distribution which is reliable for all values 
of C- We achieve this by matching to a fixed-order (FO) calculation. Although a complete 
phenomenological analysis would require matching to a next-to leading order (NLO) QCD 
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Figure 3: Matching of the NLL resummed distribution to the LO one for R = 0.6 (on the left) 
and R = 1.0 (on the right). 



calculation, for the purpose of this paper we limit ourselves to LO matching. We compute 
the differential jet-mass distribution at 0(a s ) using Comix. The result is plotted in Fig. || 
(dashed black line): the differential distribution -^§^ diverges logarithmically in the small- 
mass limit. The dotted green line instead represents the resummed result. The matched 
curve (shown in solid red) is obtained by straightforwardly adding the two contributions 
and removing the double-counted terms, i.e. the expansion of the resummation to 0(a s ): 



1 tfo~NLL+LO _ 

a <i In £ a 



dlnC dln( dlnC 



(5.1) 



We note that the matched result coincides with the resummation at small £, because 
the logarithmically divergent contributions to the LO distribution are cancelled by the 
expansion of the resummation. Moreover, the matched distribution follows the LO order 
one in the opposite limit. In particular, we note that the LO result exhibits an end-point: 



■> rn 2 j 2p t k t 2p t k t 

C = ~T = 1 —T12 ( COsh V ~ COS 0) = 2 , ;2 I n 1 1 COsh V ~ COS < P) ' 5 - 2 

Pt \P t + !k\ Pt + 2p t k t cos r 

which leads to 



C max = . / max C 2 = tan | = | + 0(R 3 ). (5.3) 

y y +4 ' <-R ^ & 

This feature is not present in the resummed distribution, because the jet does not recoil 
against the emission of the eikonal gluon, but it is restored thanks to matching. 

5.3 Numerical estimate of constant term C\ 



In this section we investigate the numerical impact of the constant C\ defined in Eq. (3Tq), 
in the case of Z+jet events, where we measure the mass of the highest pr jet. Different 
contributions build up the constant terms at 0(a s ): we have one- loop virtual corrections, 
terms that cancel the renormalisation and factorisation scale dependence of the Born pieces 
and contributions with multiple partons in final state that may, or may not, end up in the 
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same jet. One- loop virtual corrections have the same kinematical and flavour structure 
as the corresponding Born subprocess 5 and, consequently, they need to be suppressed 
by the same resummed exponent. In order to cancel infrared divergencies, one needs to 
consider real emissions in the soft and collinear limit as well. The final-state singularities 
are precisely the source of the logarithms we are resumming and these configurations can 
be mapped onto one of the Born subprocesses (<5). Initial-state collinear singularities, 
which do not give rise to any logarithms of the jet mass, must be absorbed into the parton 
densities, leaving behind terms that depend on the factorisation scale. Finally, we also need 
to consider kinematic configurations where the two final-state partons are not recombined, 
resulting into Z+2 jet events, and suppress the hardest one with the appropriate exponent. 
Thus, we should perform the calculation of the NLO cross-section that appears in Eq. ( 3.18| ) 
keeping track of the kinematics of the final-state partons, in order to separate between g 
and q components. Although this can be clearly done by computing these corrections 
analytically, for the current analysis we use the program MCFM [54], which does not 
trivially allow us to do so. Nevertheless, we are able to compute the finite part of the 
virtual corrections, for the different initial-state channels. We suppress virtual corrections, 
as well as the integrated term in Eq. (3.18) with their appropriate form factor, f( s \ We 
then multiply all the remaining real corrections by either a gluon or a quark form factor, 
producing a band. 

Our findings are plotted in Fig. ^j, for R = 0.6, on the left, and R = 1.0, on the right 4 . 
When C\ is included, we normalise the distribution to the total NLO rate, rather than the 
usual Born cross-section. In order to avoid large NLO corrections [55], we put a cut on the 
Z boson transverse momentum prz > 150 GeV. We have found that this leads to X-factor 
K = 1.45 and K = 1.57, for R = 0.6 and R = 1.0, respectively. With this set of cuts, the 
impact of C\ is not too big, but definitely not negligible. The complete calculation of this 
constant is therefore necessary in order to be able to perform accurate phenomenology. 



4 We have suppressed C\ with the full global resummed exponent, producing terms beyond our NLL 
accuracy, which we do not control. A more precise analysis would involve the complete determination of 
Ci, suppressed only with double- logarithmic terms, together with an uncertainty band, assessing the impact 
of higher logarithmic orders. 
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Figure 4: The impact of the NLL constant C\, for R — 0.6 jets (on the left) and R — 1.0 (on the 
right). The band is produced by suppressing the real radiation contributions with a quark or gluon 
jet form factor, as explained in the text. 



5.4 Comparison to Monte Carlo event generators 

In this section we compare our resummed and matched result NLL+LO to three standard 
Monte Carlo event generators: Sherpa [52], Pythia [56] and Herwig++ [57]. Monte Carlo 
parton showers are powerful tools to simulate complicated final states in particle collisions. 
When interfaced with hadronisation models, they are able to describe the transition be- 
tween partons and hadrons. Moreover, they provide events which are fully differential in the 
particles' momenta. Also they provide models of other non-perturbative effects at hadron 
colliders such as the underlying event where it is not possible to have phenomenological 
estimates from first principles of QCD. As stated before however, despite their usefulness 
and successes, it is quite difficult to assess the theoretical precision of these tools. For 
this reason, comparisons between parton showers and analytic calculations, which have a 
well-defined theoretical accuracy, form an important part of QCD phenomenology. For the 
case of jet-mass it has been noted in recent ATLAS studies [12] for jet masses that Pythia 
with hadronisation and the underlying event gives a reasonably good description of the 
data. We would therefore expect our resummation to be in some accordance with Pythia 
though we should stress that we do not include any non-perturbative effects. Hence we 
compare our results to the parton shower aspect of the various event generators on its 
own. While this is in principle correct, in practice one should be aware that there can be 
considerable interplay between the shower level and the non-perturbative effects in various 
event generators so that these programs may only return more meaningful results when all 
physical effects (perturbative and non-perturbative) are considered together. We should 
bear this caveat in mind while attempting to compare a resummed prediction with just the 
parton shower models in event generators. 

The results of the comparison are shown in Fig. H Our NLL+LO result for R = 0.6 
is shown in red (the band represents the uncertainty due to the incomplete determination 
of C\). The Monte Carlo results are obtained with the same parton densities as in the 
resummed calculation and the same set of cuts. For Sherpa (blue line) and Pythia 
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(green line) we fix [Mr = fj,p = 200 GeV, while for Herwig++ (magenta line) we use the 
default transverse mass of the Z boson. At the shower level, Sherpa and Herwig++ 
appear to perform quite similarly. They produce fairly broad distributions, which are not 
very much suppressed as £ — >• 0. The distribution obtained from Pythia instead produces 
a curve which is much closer to our resummed result. Although the position of the peaks 
differ by 5£ = 0.01 (5mj ~ 2 GeV), the height and general shapes appear in agreement. 

The agreement between the different Monte Carlo generators is restored when hadro- 
nisation corrections are included, as demonstrated in Fig. |(| Pythia and Herwig++ 
produce very similar results, while the distribution obtained with Sherpa is broader, but 
not too different. Clearly, in order to compare to collider data, one must also include the 
contribution from the underlying event. 

It is also interesting to compare the resummed prediction we computed in this paper 
with a shift to the right to account for hadronisation corrections to the event generator 
results after hadronisation. The shift approximation, initially suggested in [58], should be 
valid to the right of the Sudakov peak but will certainly break down in the vicinity of 
the peak (see also the discussion in [26] and references therein). The amount of the shift 
is related to the non-perturbative correction to mean values of jet masses derived in [10] . 
From that reference, we note that the v = m 2 j/p^ distribution should be shifted by an 
amount aR/px (the correction to the mean value) with a dimensionful coefficient a one 
can take to be of order of a few times Aqcd- The results of our calculations with a non- 
perturbative shift are shown in Fig. ||. From there one notes that a shift of 1.5 GeV -R/p*r 
(for the v distribution which we carry over to the £ distribution plotted in Fig. ^) where we 
take pt to be the value of the lower bound (200 GeV in this case) on transverse momentum, 
yields an excellent agreement with the Herwig+- |- result after hadronisation. On the other 
hand, a slightly larger shift of 2.0 GeV R/pt yields a good agreement with Pythia. We 
have truncated the shifted resummed result near the peak of the resummed distribution, 
as we would expect the shift to not be meaningful beyond this region at the very best. 
Although we have made a crude estimate of hadronisation effects, and one may be able to 
better compute these corrections, we do note that within non-perturbative uncertainties our 
results are compatible with the most widely used event generator models. We therefore 
anticipate that after the further improvements we have in mind are accomplished, our 
results may be directly used for phenomenological purposes. 
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Figure 5: Comparison of our resummed and matched result NLL+LO (in red) to standard Monte 
Carlo event generators, at the parton level. 
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Figure 6: Results for the £ distributions obtained with standard Monte Carlo parton showers, with 
hadronisation corrections (dashed lines) compared to analytical resummation with non-perturbative 
shifts (shaded bands) as explained in the main text. 



6. Dijets at the LHC 

In this section we provide numerical predictions for the jet mass distribution in dijet events. 
As before, we consider proton-proton collision at y/s = 7 TeV, with jets defined according to 
the anti-kt algorithm [27]. The main complication with respect to the Z+jet case previously 
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Figure 7: The NLL+LO jet mass distribution for dijets, with different approximation to the 
resummed exponent. 



discussed is a more complicated colour algebra, which results into a matrix structure of 
large-angle soft gluon radiation. In order to simplify our resummed calculation, we work 
at fixed kinematics, i.e. we demand the jets' transverse momentum to be pt = 200 GeV 
and their rapidity separation to be |Ay| = 2 (at Born level we only have two jets). We 
remind the reader that we consider 



Ida 

ad( 



1 f da da 
a \d(i d( 2 



(6.1) 



Ci=fe=C 



We match the resummation to a LO calculation of the jet mass distribution obtained 
with Nlojet++ [7], according to Eq (|5TT|) . In Fig. we plot our NLL+LO result with 
different approximation for the resummed exponent Eq. ( 4.1 3| ) : jet-functions (blue), with 
the inclusion of 0(R 2 ) corrections (green) and non-global logarithms (red), as explained 
in detail for the Z+jet case. Although the corrections to the jet-function approximation 
are less pronounced than in the Z+jet case, they are still sizeable and must be taken into 
account. In our understanding, the perturbative part of the resummed result of Ref [19,20] 
is precisely the one captured by the jet-functions (plus an approximated treatment of the 
one-loop constant C%). 

In order to obtain a more realistic prediction, one would need to integrate over the 
appropriate cuts in transverse momentum and rapidity and match to NLO, which, in 
principle, should not pose any difficulties. More delicate is instead the determination of 
the constant C±, although this issue has been addressed in Ref. [16,38] for the case of global 
event-shapes. 
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7. Conclusions and Outlook 

In this paper we have provided resummed expressions to NLL accuracy for jet-mass distri- 
butions in hadron-hadron collisions, for jets of arbitrary radius and defined in the anti- fa 
algorithm. In particular, we have considered Z boson production in association with a jet 
and jet masses in dijet production. We have improved upon existing studies of jet masses 
and shapes (many of which use e + e~ annihilation as a model for jets produced in hadron 
collision processes) by incorporating initial state radiation effects as well as accounting for 
non-global logarithms in the large N c limit. We have matched our results to leading-order 
predictions to account for non-logarithmic terms to this order and have commented on the 
role of the coefficient function C\a s (corrections of order a s to the basic resummation off the 
Born configuration). Finally, we have compared our results to all the leading Monte Carlo 
event generators Pythia, Sherpa and Herwig++, with and without non-perturbative 
hadronisation corrections. 

We have found that firstly ISR and non-global logarithms play an important role even 
at relatively small values of jet radius such as R = 0.6 and hence cannot be ignored while 
discussing inclusive jet shapes at hadron colliders. Our calculations for non- global loga- 
rithms both at fixed-order and at all orders represent to our knowledge the first attempt 
to calculate these terms beyond the simpler cases of e + e~ annihilation and DIS. On com- 
paring our results to event generators widely used for phenomenology we find that at the 
purely perturbative level the best agreement is with the Pythia shower, with an apparent 
small shift accounting for much of the difference between the analytical resummation and 
the parton shower estimate. The differences with Sherpa and Herwig++ on the other 
hand are more marked especially towards smaller jet masses as one approaches the peak of 
the distributions. After hadronisation corrections are applied in the event generators and 
also applying a shift to the analytical resummation to account for hadronisation, we are 
able to obtain very good agreement with Pythia and Herwig++ with slightly different 
shifts required in either case, for jet mass values to the right of the Sudakov peak of the 
distribution. For smaller jet masses we do not expect a simple shift of the analytical re- 
summation to reproduce the correct result and we observe here discrepancy with all Monte 
Carlo generators which is to be expected. However, this region of very small jet masses 
of the order of a few GeV will not be of interest in LHC phenomenology in any case. We 
may thus expect that our results should be of direct phenomenological value even pending 
some of the improvements we intend to make in the near future. 

For the immediate future we aim to improve our results chiefly by taking proper ac- 
count of the order a s coefficient function C\ and computing the various pieces of C\ which 
originate from different regions and suppressing these by an appropriate form factor rather 
than the cruder treatment reported in the text, which was aimed at producing an uncer- 
tainty band associated to a lack of correct treatment of C\ . When this is done we will be in 
a position to carry out an NLO matching and estimate the uncertainty of our theoretical 
calculations accurately, which will be important in the context of phenomenology. 

We have not taken directly into account non-perturbative effects in this paper, incor- 
porating them as a simple shift of the perturbative spectra as for the case of global event 
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shapes [37]. We can study non-perturbative corrections in more detail using the methods 
outlined in [10], but for the moment we note that the Monte Carlo event generators we 
have studied, contain differences in their estimate of hadronisation which should be ex- 
plored further along with studying the underlying event (UE) contributions. Since our 
predictions are valid for any value of jet radius R, one can hope that phenomenological 
studies selecting a small value of R would help in better isolating the perturbative con- 
tributions, since the hadronisation and UE corrections for rrVj/p^ scale respectively as R 
and R 4 [10]. On the other hand moving to a larger R after pinning down the perturbative 
content would help to more accurately constrain the non-perturbative models in various 
approaches. 

Additionally, although in this article we have treated a single variable, the jet mass, it 
is actually straightforward to extend our treatment to the entire range of angularities as for 
instance were explored for jets in [17, 18]. The basic calculations we carried out here can 
be easily extended to include those variables as well as any variants of the jet mass itself 
such as studying the jet mass with an additional jet veto. Once we have accomplished an 
NLO matching we therefore intend to generalise our approach to accommodate a range of 
substructure variables in different hard processes. We hope that our work will eventually 
lead to improved estimates of the accuracy and hence more confidence in a detailed under- 
standing of jet substructure than is the case presently, which could in turn be important 
for a variety of substructure applications at the LHC. 
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A. Dipole calculations for the global part 

In order to carry out the calculations for the various dipole terms required for the resum- 
mation detailed in the main text we shall make use of the following Born kinematics for 
the hard incoming and outgoing particles, 

pi = ^X! (1,0, 0,1), 

P2 = ^x 2 (1,0, 0,-1), 
P3 = Pt (cosh y, 1, 0, sinh y) , 

P4 = Pt (cosh y r , -1, 0, sinh y r ) , (A.l) 
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where, as in the main, text p\ and p 2 denote the four- moment a of the incoming hard 
partons, ps that of the jet whose mass is being considered and p$ the four-momentum 
of a recoiling jet. In the case of a final state vector boson instead of a recoiling jet the 
contribution from p^ to the dipole calculations is of course absent as the Z boson is colour 
neutral. Further we have indicated by x\ and X2 the fractional momenta of incoming 
partons relative to the colliding hadrons Pi and P2 respectively. Also y denotes the rapidity 
of the measured jet, y r that of the recoiling jet and pt their transverse momenta. One can, 
by ignoring recoil, use the Born kinematics outlined above even in the presence of a soft 
emission, whose momentum can be expressed as 

k = kt (cosh r], cos (ft, sin 0, sinh 77) ; (A. 2) 

where transverse momentum kt, rapidity r\ and azimuth cp are all measured with respect 
to the beam direction. 

We shall now set up the various dipole contributions to the jet mass cross-section. 
We recall from the main text that eventually these dipole contributions will exponentiate 
with appropriate colour factors, so what follows below is in essence a calculation of various 
components of the resummed exponent. 

A.l In-in dipole 

The in-in dipole is made from incoming partons 1 and 2 and provides the contribution 

K 12 = J k t dk tdr} - ^ {pik)[p ^f [~ (coahfo - y) - cos0) - vj 

e( J R 2 -((r / -y) 2 + 2 )) (A.3) 

We next note that k 2 12 = 2 ^ Pl ^^ ) ,fc - ) = k 2 and carry out the integration over rj and 
4>. In doing so we can neglect the dependence of the jet mass on 77, <p and just retain its 
dependence on kt since the coefficient of kt which contains the rj, (j) dependence, produces 
terms that are below single logarithmic accuracy. Performing the integral over r\ and 4> 
over the interior of the jet region gives us to the relevant single- log accuracy that 



= p 2 r ^m^, (a.4) 



K12 

Jvpt 

where the lower limit of integration stems from the constraint on the jet-mass variable. 
We note that the dipole consisting of the two incoming partons gives rise to a pure single- 
logarithmic behaviour as is evident from the fact that the emitted gluon is inside the jet 
region, away from the hard legs constituting the dipole and hence there are no collinear 
enhancements. The soft wide-angle single logarithm we obtain is accompanied by an R 2 
dependence on jet radius, reflecting the integration over the jet interior. 

A. 2 In-recoil dipoles 

As before for the in-in case, for a dipole made up of an incoming parton and one recoiling 
against a measured jet, we have that the emitted gluon is in the interior of the jet away 
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from the emitting legs, and produces no collinear logarithm. We have that 

<r> I i at, a d( f )as (Pi-Pi) r^f 2kt ( \ A\ 

TZ U = / hdktdri- — —9 — (coshfa - y) - cos <p) - 

J 2vr 2vr {Pi-k){p4.k) \ p t 

G( J R 2 -((r ? -y) 2 + 2 )). (A.5) 

However , in this case we have that 

«J 14 = 2{pi - k){Pi - k) = 2Af (coshfr - y r ) + cos 0) . (A.6) 
To single logarithmic accuracy, we can write the result into the form 

/* -a s (xp t ) [ e ""* 6 (R 2 - ((v - V? + 2 )) • (A.7) 

J v x J cosh yq — y r ) + cos <p 

It is possible to evaluate the integral over 7], <fi as a power series in R, such that the overall 
result for this dipole reads 

fUx asjxpt) (R 2 e A y R 4 h 4&y^ n( ^\ 

L VTa + coshA,) + 32 sech t + ° {R } ) (A - 8) 

where Ay = y — y r For the 7^-24 dipole one obtains the same result as above, with the 
replacement of Ay by — Ay. 

A. 3 Jet-recoil dipole 

For this dipole we have also an involvement of the parton that initiates the measured jet 
as one of the hard emitting legs. The consequent collinear singularity, in addition to the 
pole due to soft emission, gives rise to double-logarithmic contributions. In this case one 
needs to evaluate the argument of the running coupling more carefully whereas for soft 
wide-angle pieces, dealt with previously, all that is important at our accuracy is that the 
argument of the running coupling is proportional to the energy of the soft gluon being 
considered. 

In the present case we start with the integral 

7<-34 = / k t dk t dr]- — —9 — (cosh^ - y) - cos 0) - v) 

J 2vr 2vr {p3-k)(pA.k) \ p t ) 

G( J R 2 -((r ? -y) 2 + 2 )) (A.9) 

where 

2 _ 2(p 3 .k)(p 4 .k) _ 2 (cosh (y r - rj) + cos <t>) (cosh (y - if) + cos 4>)) 
K t,34 — 7 \ — iK t 1 : m • i,A.iu; 



Next we introduce the variables r and 8 such that y — r/ = r cos 9 and <f> = r sin 9 such 
that the in jet condition (y — r/) 2 + eft 2 < R 2 , simply produces r < R, while 9 takes 
values between and 2n. Thinking of Eq. A.9 as a power series in r we can isolate the 
collinear-enhanced terms which produces double logarithms and separate it from terms 
that are finite as r — > 0, which upon integration produce terms that behave as powers of R 
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corresponding to soft wide-angle contributions. The contribution 7^34 can thus be broken 
into terms that contain collinear enhancements (and hence leading double logarithms) and 
pure soft single- logarithmic pieces: 

n u = nt+nt, (a. 11) 

where the pure soft terms can be expressed as a power series in R: 

,soft _ f 1 us(xpt) dx ( R 2 2 Ay R 4 ,4 Ay , . , , A , a , ^ 



nf " = J v 2^~{T taUh 2 + 1152 Sech 2 (- 5 + coshA ^ +<W 

(A.12) 

where we have neglected terms that vanish as v — > 0. 

The collinear contributions require us to be more precise in treating the argument of 
the running coupling. In particular, in the limit of small r, which is the collinear region, 
we obtain 

^mii f R dr 2 d9 I' 1 dx a s , . . „„, 

^34 1L =/ — T- / — o 1 Osm) A.13 
J v r 2 2vr Jv x 2ir 

where we have written the result in terms of a dimensionless energy fraction x such that in 
the collinear limit r <C 1, one has that nt,zil { r Pt) = x - The lower limits of the energy and 
angular integrals follow from the restriction on the jet mass to be greater than v, while the 
upper limit of the angular integral over r 2 corresponds to the jet boundary. 

It is straightforward to obtain the result corresponding to the fixed-coupling limit from 
Eq. flPaD , which gives: 



The running coupling result can be compared directly to that obtained in the e + e~ 
case in our previous work [23]. In order to do this we introduce the transverse momentum 
with respect to the emitting jet kt j = xrpt, and changing variable we obtain 



t,J / 
(A.15) 

Carrying out the integral over r, we are left to evaluate the integral over kt j which 
reads 



1-> coll. 
/-c 34 




V 

R 2 



(A.16) 
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A. 4 In-jet dipoles 

Now we consider the dipole formed by the triggered jet and one of the incoming partons. 
Since here too the measured jet is involved as one of the legs of the hard dipoles, we can 
carry out precisely the same manipulations as for the IZ34 dipole immediately above. Doing 
so yields the same soft-collinear piece as before 

nif- =n c 2 f=nit, (A.i7) 

while the soft wide-angle piece reads 

= mf = jf 1 «=£p> ^ (t + W8 + ° ^) • (A - 18) 

B. Dipole calculations for the non-global contribution 

We now turn our attention to the evaluation of the different contributions arising from 
correlated-gluon emissions. 

B.l Dipoles involving the measured jet 

We start by considering the the three dipoles involving the triggered jet I13, I23 and 
The best way to evaluate them is to boost into the frame where the two legs of the dipole 
are back to back: 

p 3 = §(1,0,0,1), (B.l) 

Pi = |(1, 0,0,-1), i = 1,2,4, (B.2) 

where Q is a function of the components of the original momenta and pi (Q 2 = 2pz.pi). 
In this frame the legs of the dipoles are sitting at infinite rapidities, so we have to write 
the jet-algorithm 6 function in terms of the polar angle with respect to the jet. We have 



'i3 



-fdrn^f d m f. - ( cos ^-^) e (e R - e 2 ) e (* - e R ) , (b.s) 

J 2tt J 2ir cosh(??i - 772) - cos(0i - 2 ) 



where Or denotes the jet's opening angle. As previously noticed in the literature, we have 
the result is the same as in the hemisphere case: 



/oo 1 — In tan ^ 2 7T ' 2 



B.2 The remaining dipoles 

Let us consider the dipole from 2 incoming legs I\ 2 . 
We have to evaluate 



T A f A d ^ f A d< ^ 2 COs(0l-</> 2 ) , 2 . 2 2 , 

h2 = 4 / dru— / df] 2 - T7 ; u —t® (R - (r? 2 + <f> 2 )) 

J 2ir J 2tt cosh(r/i - r? 2 ) - cos(0i - <p 2 ) v ' 

9 {rjl + 4>\- R 2 ) . (B.5) 
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Let us change the variables to 

V1-V2 = A cos (3 , 0i - 02 = A sin B 

772 = pcosa, 02 = psina. (B.6) 

Doing so and taking care of the fact that — ir < <f>i < ir we get 



/ 



6 [-7T — (psina + Asin/3)] [-7T + (psina + Asin/3)] pdp—XdX— (B.7) 

27T 27T 

The step function implies 

r>2 _ 2 _ \2 

cos(a - /3) > (B - 8) 

This is automatically satisfied if (R 2 — p 2 — X 2 )/(2pX) < —1. It is never satisfied if (R 2 — 
p 2 — A 2 )/(2pA) > 1. Hence, the step function in question is only active when — 1 < 
(R 2 — p 2 — A 2 )/(2pA) < 1. For a non-zero result we must then have p + A > R. 

Let us consider first the region of integration where the first step function is active. 
We have two integrals 

f R , ,> cos(Asin/3) dB f R , da ^ , 9 9 ^ N 9 . 

/ AtZA — — ± ^TT-^^r / Pdp-z-® {p + A + 2pA cos (a — 13) — R ) 

J cosh(Acosp) - cos (A sin p) 2tt Jr_\ 2ir v ' 

f 2R , „ cos(Asin/3) dB f R , da ^ , 2 >2 , „ 2X 

+ / AdA — — ^ — — / pdp — G(p 2 + A 2 + 2pAcos a-/3)-i? 2 ), 

7 R cosh(Acos/3 - cos As in/3 2vr 2vr VP K v ' ' ; ' 



(B.9) 

where we have not written the other step functions that are automatically satisifed since A 
and p are small here. The first integral yields « 0.61.R 2 while the second one is « 0.37-R 2 . 

When the cosine condition is automatically satisfied we have the following range of 
values 

R<X<2R,0<p<X-R (B.10) 

and 

A > 2R, < p < R. (B.ll) 

The first integral is 

f 2R , „ cos(Asin/3) d/3 I"*-* , da 

/ AdA — — ^ ^7- — — / pdp — (B.12) 

J R cosh(A cos B) — cos(A sinp) 2ir Jq 2tt 

Note that for small R and the above range of values of p all other theta functions are also 
satisfied as e.g it is always true that ir > psina + Asin/3. One can evaluate the integral 
above easily and the result is « 0.19i? 2 . 
The remaining integral is 



/ 

J 2 



00 , „ cos(Asin/3) dB f R , da 

XdX — — — — — — / pdp — 

2 r cosh(Acosp) — cos (A sin pj 2tr J 2tt 

6 [7r - (psina + Asin/3)] 9 [ir + (psina + Asin/3)] (B.13) 
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We note that we can split the A integral to go from 2R to unity and then unity to infinity. 
The integral from 2R to unity produces the following result: 

-R 2 ln2i?-0.12i? 2 

which leaves to evaluate the integral 



„ cos(Asin/3) d(3 f R , da 

XaX — — — — : — — — / pdp- 



1 cosh(Acos/3) — cos(Asin/3) 2ir J 2ir 

6 [vr - (psina + Asin/3)] 9 [tt + (psina + Asin/3)] (B.14) 

This goes as kR 2 with k a constant to be determined. To evaluate k it proves simplest to 
first take the derivative with respect toR 2 which turns the theta function involving p into 
a delta function. The p integration is then trivial and the rest of the integral can be done 
numerically and yields k = 0.12. We thus have that the overall result is 0.12i? 2 . Note that 
it is interesting that this cancels the — 0.12R 2 obtained previously. Hence the full result is 

/12 ~ 4 [1.17R 2 - R 2 ln2R-0.12R 2 + 0.12R 2 } +0(R 4 ) « 4 [l.l7i? 2 - i? 2 ln2ii] +0(R 4 ). 

(B.15) 

The above expansion well describes the full answer for I12 up to values of R ~ 0.6. For 
larger values of the jet radius we use a full numerical evaluation of the coefficient I±2- 

Next we consider the in-recoil dipole /14. The result can be expressed in terms of the 
dipole I12 and a rapidity dependent function n(Ay), which we evaluate numerically 

(l + eZV> „ (1.17R 2 - R 2 \n2R) + (1 + ^ K(Ay)R 2 . (B.16) 
(l + coshAy) 2V ; 1 + coshAy y y > y 1 

Finally, the dipole hi is easily obtained as 

I 2A = I 14 (-Ay). (B.17) 

C. Resummation formulae 

In this section we collect the explicit expressions of the function which build up the 
resummed results: 

/l(A) = -^^[a-^wi^A^u-A^na-A)], 

/2(A) = -^[21n(l-A)-ln(l-2A)] 



4vr 2 /3 

§1 
2vr/3 3 



In (1 - 2A) - 2 In (1 - A) + - In 2 (1 - 2A) - In 2 (1 - A) 



WA) = ~ln(l-A), 
47rpo 

WA) = ~ln(l-A), 
/i.a.(A) = -^ln(l-2A). (C.l) 
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A = (3oa s L, L = In ^ and a s = a s {ptR) is the MS strong coupling, which at aimed 
accuracy we need to consider at two- loops: 



a s (k t ) - 



1-p 



2 ' { 2 ,/3iln(l-p) 
1 - a s (fi )■ 



Po 1-P 



,9 = a s (^ 2 )/3 ln 



where the coefficients of the QCD /^-function are defined as 

. llC A -2n f . 17^-5(7^-3^71/ 
Po = 77; , Pi = 



12vr 



24vr 2 



and the constant K is given by 



(C.2) 



(C.3) 



K = C A 



{67 



- 



7T _ 



V 18 6 



-n 



(C.4) 



Finally, we report the coefficients a, 6, c obtained by fitting the functional form 



/«(*) 



1 + (M) 2 ^2 

1 + •• 



(C.5) 



to the numerical dipole evolution for non-global logarithms. When we consider a dipole 
involving the measured jet we find 



a = 0.85C A , b = 0.86C A , c = 1.33, (C.6) 

independently of the jet radius R. The results for the other dipoles are reported in Tables [j] 
and |2[ Note that the results for I±4 and I24 depend on the rapidity separation between the 
jets and they have been computed for |Ay| = 2. 



R 


I12 


a 


b 


c 


0.1 


0.11 


2.67 Ca 


0.93 Ca 


1.66 


0.2 


0.34 


2.27 C a 


1.08 Ca 


1.66 


0.3 


0.62 


2.05 C A 


1.14 Ca 


1.66 


0.4 


0.92 


1.99 C A 


1.31 Ca 


1.66 


0.5 


1.22 


0.97 Ca 


0.00 Ca 


1.33 


0.6 


1.52 


0.85 Ca 


0.00 Ca 


1.33 


0.7 


1.80 


0.96 Ca 


0.29 Ca 


1.33 


0.8 


2.05 


1.12 Ca 


0.56 Ca 


1.33 


0.9 


2.28 


1.12 Ca 


0.63 Ca 


1.33 


1.0 


2.49 


1.24 Ca 


0.80 Ca 


1.33 


1.1 


2.67 


1.33 Ca 


0.98 Ca 


1.33 


1.2 


2.82 


1.45 Ca 


1.15 Ca 


1.33 



Table 1: Numerical results for the coefficients that parametrize the resummation of non- global 
logarithms. 
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H 


IlA 


a 


Z, 



c 


j 


a 


1. 




c 


0.1 


0.09 


2.74 La 


0.86 Ca 


1.66 


0.00 


3.85 Ca 


r\ r\r\ 

0.00 Ca 


2.00 


0.2 


0.28 


2.57 C a 


1.27 Ca 


1.66 


0.01 


3.28 Ca 


0.31 Ca 


2.00 


0.3 


0.51 


2.27 Ca 


1.31 Ca 


1.66 


0.02 


2.84 Ca 


0.23 Ca 


2.00 


0.4 


0.77 


2.02 Ca 


1.25 Ca 


1.66 


0.03 


2.89 Ca 


0.59 Ca 


1.66 


0.5 


1.04 


1.03 Ca 


0.00 Ca 


1.33 


0.05 


2.90 Ca 


0.77 Ca 


1.66 


0.6 


1.31 


0.94 Ca 


0.00 Ca 


1.33 


0.07 


2.76 Ca 


0.81 Ca 


1.66 


U. ( 


1.0 < 


u.yo ua 


U.14 Oa 


i.oo 


U.1U 


Z.OO (_/A 


U. /O OA 


l.DD 


0.8 


1.81 


0.96 C a 


0.31 Ca 


1.33 


0.13 


2.45 Ca 


0.78 Ca 


1.66 


0.9 


2.04 


0.86 Ca 


0.21 Ca 


1.33 


0.16 


2.50 Ca 


0.93 Ca 


1.66 


1.0 


2.25 


0.90 Ca 


0.30 Ca 


1.33 


0.20 


2.45 Ca 


0.99 Ca 


1.66 


1.1 


2.44 


1.18 Ca 


0.72 Ca 


1.33 


0.24 


2.39 Ca 


1.05 Ca 


1.66 


1.2 


2.61 


1.00 Ca 


0.51 Ca 


1.33 


0.30 


2.56 Ca 


1.30 Ca 


1.66 



Table 2: Numerical results for 
logarithms. Note that the above 



the coefficients that parametrize the resummation of non-global 
results are valid in the case |Ay| = 2. 
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